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THE CHEMICAL SIGNATURE OF A RELIC STAR CLUSTER IN 
THE SEXTANS DWARF SPHEROIDAL GALAXY — IMPLICATIONS FOR NEAR-FIELD COSMOLOGY 
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ABSTRACT 

We present tentative evidence for the existence of a dissolved star cluster at [Fe/H] = —2.7 in the Sextans 
dwarf spheroidal galaxy. We use the technique of chemical tagging to identify stars that are highly clustered 
in a multi-dimensional chemical abundance space (C-space). In a sample of six stars, three, possibly four, 
stars are identified as potential cluster stars. The initial stellar mass of the parent cluster is estimated from 
two independent observations to M* i n i t = 1.9^o'g (l.6t-Q's) x 10 5 -M©' assuming a Salpeter (Rroupa) initial 
mass function (IMF). If corroborated by follow-up spectroscopy, this star cluster is the most metal-poor system 
identified to date. Chemical signatures of remnant clusters in dwarf galaxies like Sextans provide us with a very 
powerful probe to the high-redshift Universe. From available observational data, we argue that the average star 
cluster mass in the majority of the newly discovered ultra-faint dwarf galaxies was notably lower than it is in the 
Galaxy today and possibly lower than in the more luminous, classical dwarf spheroidal galaxies. Furthermore, 
the mean cumulative metallicity function of the dwarf spheroidals falls below that of the ultra-faints, which 
increases with increasing metallicity as predicted from our stochastic chemical evolution model. These two 
findings, together with a possible difference in the ([Mg/Fe]) ratio suggest that the ultra-faint dwarf galaxy 
population, or a significant fraction thereof, and the dwarf spheroidal population, were formed in different 
environments and would thus be distinct in origin. 

Subject headings: Galaxies: dwarf- Galaxies: formation -Galaxies: star clusters: general - Galaxy: formation 
- Stars: abundances - Stars: Population II 



1. INTRODUCTION 

The formation and long-term evolution of dwarf galaxies are 
topics of great interest today. A substantial fraction of the 
dwarf galaxy population of the Local Group is expected to 
have already merged with the halos of M3 1 and the Galaxy. 
This is supported by the evidence that the stellar Galactic Halo 
has several distinct chemical similarities that it shares with 
the existing population of dwarf galaxies (Kirby et al. 2008; 
Tolstoy et al. 2009 and references therein; Starkenburg et al. 
2010; Frebel et al. 2010a, b; Norris et al. 2010a, c). Despite 
these similarities, the systems are not chemically identical and 
exactly how the merging process progressed, how the dwarf 
galaxies that survived are connected to their disrupted coun- 
terparts, and how these still-surviving galaxies relate to each 
other in a cosmological context, is not fully understood. 

The study of old stellar populations in these galaxies 
provides us with an indirect link to the high-redshift Uni- 
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verse. All Local Group dwarf galaxies, including the re- 
cently discovered ultra-faint dwarfs, are known to contain an- 
cient, metal-poor stars (Grebel & Gallagher 2004; Tolstoy et 
al. 2009). The star formation efficiency in dwarfs is very 
low which means that there have been relatively few star- 
formation events over cosmic time. The effects of stochastic- 
ity therefore plays a relatively big role (Audouze & Silk 1995; 
Karlsson 2005; Karlsson & Gustafsson 2005). This aspect 
must be addressed when seeking to understand the evolution 
of dwarf galaxies, particularly at early times. 

A key characteristic of star formation is the creation of 
star clusters. The initial cluster mass function (ICMF), if it 
can be determined reliably, provides important information 
on the physical conditions under which the clusters formed 
(Elmegreen 2010; Larsen 2009). A flat ICMF appears to 
be indicative of high-pressure environments (e.g. starbursts, 
mergers), whereas a steep ICMF is typically associated with 
quiescent environments (e.g. solar neighbourhood). At the 
present time, the most metal-poor star clusters known have an 
iron abundance just below [Fc/H] = —2. One of the globu- 
lar clusters (Cluster 1) in the Fornax dSph currently holds the 
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Figure 1. Chemical abundance ratios of six very metal-poor stars in the Sextans dSph. Abundance data are taken from Aoki et al. (2009). Stars are color coded 
according to their location in the [Mg/Fe] - [Fe/H] diagram. The three stars with nearly identical [Mg/Fe] (colored in different shades of blue) also have very 
similar [Fe/H]. These stars all clump together in Ti, Cr, and Ba, as well. One of the "blue" stars, S 10 — 14 (dark blue), is substantially deficient in Na. The 
scandium abundance could not be measured in S 10 — 14, while nickel could not be measured in S 10 — 14 and S 14 — 98. In the [Mg/Fe] (upper, middle left) 
and [Ti/Fe] (lower, outer left) diagrams, the medium (S 14 — 98) and light blue (S 11 — 13) markers, respectively, are made smaller for visual purposes. 



record with a metallicity [Fe/H] = -2.5 (Letarte et al. 2006). 
But not all clusters have survived as gravitationally bound ob- 
jects to the present epoch. Much like the halo, dwarf galaxies 
have stars with metallicities well below [Fe/H] = —3 (Kirby 
et al. 2008; Starkenburg et al. 2010). The relatively simple 
environments of the low mass dwarfs raises the prospect of 
identifying disrupted star clusters at much lower metallicity 
than has been possible before. By extension, this also gives 
us a unique tool to probe the formation and early evolution 
of the present dwarf galaxy population and to determine the 
sequence of events that led to the formation of the stellar Halo 
of the Milky Way. 

To date, all open clusters have been shown to be extremely 
homogeneous in elements heavier than Na (e.g. De Silva et 
al. 2006, 2007). The same holds true for most globular clus- 
ters, particularly in [Fe/H] and heavier elements (Sneden et 
al. 2008). Bland-Hawthorn et al. (2010b) provide a theoret- 
ical framework for determining the mass of a homogeneous 
star cluster. Clusters with central column densities ~ 0.3 g 
cm~ 2 (e.g. open clusters) are homogeneous up to 10 5 M.q. 
More centrally concentrated clusters, with column densities 
~ 3 g cm -2 (e.g. globular clusters), are homogeneous up 
to at least 10 7 7W©. We can look for signs of clumping in 
chemical abundance space (C-space), particularly in heavier 
elements, in order to identify systems that have long since dis- 
solved (Bland-Hawthorn et al. 2010a). Globular clusters are 
distinguished by more complex signatures in lighter elements, 
in particular, an anticorrelation between O/Mg and Na/Al. 

In §2, we revisit the recent observations of Sextans pre- 
sented by Aoki et al. (2009). We reconsider their excellent 
analysis and demonstrate the possible existence of a disrupted 
stellar cluster at [Fe/H] = —2.7, well below the metallicity 
of the most metal poor globular cluster to date. In §3, we de- 
termine the statistical significance of the cluster in the context 
of inhomogeneous chemical evolution models. In §4, we es- 
timate the mass and discuss the likely nature of this ancient 
star cluster before discussing implications for near-field cos- 
mology in §5. Some final remarks are given in fj6] 



2. OBSERVATIONAL EVIDENCES 

In the following section, we will discuss three key observa- 
tions suggesting that a significant fraction of the most metal- 
poor stars in the Sextans dSph were once members of a mas- 
sive star cluster. 

2.1. Chemical abundance ratios 

Based on high-resolution (71 ~ 40, 000) spectroscopy, Aoki 
et al. (2009) determined the chemical abundances of six very 
metal-poor stars ([Fe/H] < —2.5) in the Sextans dSph. Four 
of these stars display a subzero [Mg/Fe] ratio, with a small 
scatter around the weighted mean ([Mg/Fe]) = —0.06 (see 
Fig. [TJ. Thus, they do not show the Mg-to-Fe enhancement 
(relative to the Sun) commonly observed in Galactic halo stars 
of similar metallicity (Cayrel et al. 2004; Barklem et al. 2005; 
Lai et al. 2008). This group of stars have very similar metal- 
licity ([Fe/H] ~ —2.7) and a closer inspection of Fig. Q]shows 
that they more or less tightly clump together in Cr and Ba as 
well (see also Fig. [2]). The clumping in the [Ba/Fe]-[Fe/H] 
plane around [Ba/Fe] ~ —1.3 (excluding S 11 — 37, see be- 
low) is particularly striking. 

At least three of the four stars also group together in Na, 
Ca, and Ti (i.e., Ti II). The star S 14 - 98 has a higher Ca 
abundance as compared to the other three subzero [Mg/Fe] 
stars, but as Aoki et al. (2009) points out, this may simply be 
due to a larger observational uncertainty in Ca for this star. 

In terms of the chemistry, the star S 11 — 37 (color coded 
green in Fig. [D may be regarded as a borderline case, and it is 
unclear whether it actually should be affiliated with the same 
group as S 10 - 14, S 11 - 13, and S 14 - 98 (color coded 
blue). It is more than 2cr away from the weighted mean of the 
"blue" stars in Ti, ler away in Ba, and slightly more than la 
away from the two "blue" stars for which Sc abundance could 
be measured (Sc was not measured in S 10 — 14). In all cases 
(including Mg), the star S 11 — 37 has an excess abundance 
as compared to the "blue" stars, and is closer to the typical 
abundances found in Galactic halo stars (e.g. Cayrel et al. 
2004). Such a systematic behavior would not be expected if 
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Figure 2. The data of Aoki et al. (2009) for Mg, Cr, and Ba displayed over the background of metal-poor stars of a number of dSphs (dark-grey symbols) and 
ultra-faint dwarf galaxies (light-grey symbols). The dSphs are: Sextans (circles, Shetrone et al. 2001; Tafelmeyer et al. 2010), Sculptor (squares, Shetrone et al. 
2003; Geisler et al. 2005; Frebel et al. 2010a; Tafelmeyer et al. 2010), Ursa Minor (diamonds, Shetrone et al. 2001 ; Sadakane et al. 2004; Cohen & Huang 2010), 
Draco (up-faced triangles, Shetrone et al. 2001; Fulbright et al. 2004; Cohen & Huang 2009), Fornax (right-faced triangles, Letarte et al. 2010; Tafelmeyer et 
al. 2010), and Carina (pentagons, Shetrone et al. 2003; Koch et al. 2008). The ultra-faints are: Ursa Major II (squares, Frebel et al. 2010b), Coma Berenices 
(diamonds, Frebel et al. 2010b), Bootes I (up-faced triangle, Norris et al. 2010c), Leo IV (right-faced triangle, Simon et al. 2010), and Segue 1 (pentagon, Norris 
et al. 2010a). Note the evident clumping of the "blue" Sextans stars in all three diagrams. Note also their anomalous positions, particularly in Mg and to some 
extent in Cr. To be consistent with Aoki et al. (2009), the background data are converted to the solar abundance scale of Asplund et al. (2005). 



the abundance differences were due to random errors alone. 
In what follows, we assume that S 11 — 37 does not belong to 
the "blue" group of stars. However, higher S/N spectra and 
a detailed investigation of the abundance response to stellar 
atmosphere parameters are necessary in order to determine 
the true origin of this star. 

It is clear from Figs. [T] and |2] that the "blue" stars tend to 
clump together in C-space. This clumping implies a common 
chemical history, and we will here postulate that these stars 
were once members of a chemically homogeneous star cluster 
that is now, at least partly, dissolved. Note that the "blue" stars 
could in principle have been distributed anywhere along the 
mean trends as defined by the data of a number of dSphs in the 
Local Group (Fig. |2] dark-grey symbols). The fact that they 
instead clump together at an anomalous position in C-space 
(particularly in [Mg/Fe], but also, e.g., in [Cr/Fe], see Fig. [2]i 
argues for a cluster membership in the past. A possible ar- 
gument against a common origin is that S 10 — 14 has an ob- 
served [Na/Fe] ratio nearly 3(7 below the other two potential 
cluster stars. This could, however, be a manifestation of the 
Na-O anti-correlation observed i n Ga lactic globular clusters. 
This will be further discussed in £14.21 In principle, a compre- 
hensive group-finding analysis (see Sharma & Johnston 2009 
for a presentation of the algorithm EnLink) would provide us 
with the likelihood of this being a true homogeneous clump 
in C-space (Bland-Hawthorn et al. 2010a). Unfortunately, the 



number of stars in the current sample, particularly the number 
of stars not belonging to the potential cluster, is too small and 
more data below [Fe/H] = —2.5 are required in order to give 
a reliable estimate of this probability. 

Finally, we note that the star S 12 — 28 (color coded yellow) 
has similar abundances to the "blue" stars in Ti and Cr. How- 
ever, the fact that it is over 5cr away in Ba makes the affiliation 
with the potential cluster stars improbable. 

2.2. Metallicity distribution function 

Starkenburg et al. (2010) performed a re-calibration of the 
Ca II triplet (CaT) lines and determined new metallicities 
for metal-poor stars in four classical dSphs, including Sex- 
tans, that are a part of the Dwarf galaxies Abundances and 
Radial-velocities Team (DART) program. The new data re- 
veal a "bump" in the metallicity distribution function (MDF) 
of Sextans around the inferred [Fe/H] ~ —2.9. By apply- 
ing to the DART data the presumed Ca-to-Fe enhance ment o f 
[Ca/Fe] = +0.25 (Starkenburg et al. 2010; see also gXT]), 
the "bump" appears around [Ca/H] ~ —2.65 (see Fig. [3]). 
Interestingly, this is close to the weighted mean abundance of 
Ca for the three potential cluster ("blue") stars identified in 
the data of Aoki et al. (2009). If the "bump" indicates a real 
excess of stars in the MDF of Sextans, this excess could be 
due to the presence of a relic star cluster. 
A "bump" at [Fe/H] ~ —3 exists also in the medium- 
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resolution (1Z ~ 7, 000) data of Sextans^ presented by Kirby 
et al. (201 la). They noticed that similar "bumps" are present 
in the MDFs of the Draco and Ursa Minor dSphs as well. 
They speculated that these "bumps" could be due to an early 
burst of star formation in these galaxies, but did not consider 
the particular possibility that the "bumps" may be indirect 
signs of star clusters. In fact, if stars were formed in clus- 
ters already at the earliest epochs of galaxy formation, irreg- 
ular and "bumpy" MDFs are likely to be expected, i.e., de- 
pending on the functional form of the local ICMF (see fj5). 
The irregularities should appear preferentially in regions of 
small number statistics, such as in the metal-poor tails of the 
MDFs of the least massive galaxies. If so, conclusions from 
a straightforward comparison between MDFs of individual 
dwarf galaxies should be drawn with some caution since any 
recognized difference could simply be due to the effects of 
stochasticity (see discussion in $5}. 

2.3. Kinematical substructure ? 

There is evidence for a kinematically cold substructure in the 
central region of Sextans, possibly originating from a remnant 
star cluster (Kleyna et al. 2004; Battaglia et al. 2011; see, 
however, Walker et al. 2006). Battaglia et al. (2011) identify 
a group of stars with cold kinematics at a projected distance 
of R < 0°.22 from the optical center. This substructure has 
an observed mean metallicity of [Fe/H] ~ —2.6 (Battaglia et 
al. 2011), w hich is close to that of the "blue" group of stars 
discussed in il2.11 and a common origin is not implausible. It 
should, however, be noted that the potential cluster stars iden- 
tified in the sample of Aoki et al. (2009) are located around 2 
core radii, at a projected distance of 0°.4 < R < 0°.6, well 
outside R = 0°.22. Evidently, if the "blue" stars and the stars 
associated with the central substructure do originate from the 
same disrupted (or disrupting) cluster, its former member stars 
are scattered over a significant fraction of the entire galaxy 
(Penarrubia et al. 2009). In this framework, it is close at hand 
to assume that the "bump" of stars in the MDF, in fact, is a 
manifestation of the very same cluster. In what follows, we 
will assume that this is the case and from the relative size of 
the "bump", we will estimate the total mass of stars that were 
once members of the cluster. 

We note that Sextans is unusual in having the lowest- 
density dark-matter halo of any of the classical dSph systems 
(e.g., Gilmore et al. 2007). The stars in the dSphs are likely 
to be located well within a constant density dark matter core, 
which implies that tidal stresses generally are low. For Sex- 
tans, they would be especially low because of the low density. 
This would inhibit the disruption of clusters, leading to longer 
disruption times (cf. Penarrubia et al. 2009; see also Walker & 
Penarrubia 201 1). N-body simulations of the disruption pro- 
cess in an environment typical for dSphs like Sextans should 
be performed in order to shed light on this issue. 

3. MODELING OF SEXTANS 

To illustrate what effect the presence of a massive star cluster 
could have on how C-space is populated with stars, and to 
estimate the mass of the cluster, we will model the chemical 
evolution of a small dwarf galaxy, comparable in mass to the 
Sextans dSph. 

2 Kirby et al. (2011a) make use of a particular metallicity extraction 
technique that determines an "overall metallicity" of spectrally low-resolved 
stars, rather than [Fe/H]. 



3.1. Stochastic chemical evolution model 

The model is based on the assumption that all stars originally 
are formed as members of chemically homogeneous star clus- 
ters. The mass distribution of newly formed clusters is known 
as the initial cluster mass function (ICMF), and is assumed to 
have the form 

diV/dM* = X (M*) oc M~\ (1) 

where M* denotes the cluster mass. Although the particu- 
lar situation for dwarf galaxies is not clear, especially at high 
redshifts, observations of the local Universe support a generic 
slope of 7 ~ 2 in most environments (e.g. Lada & Lada 2003; 
Fall et al. 2005, 2009; Elmegreen 2010), and cluster masses 
typically within the range 50 < M^/Mq <2x 10 5 (Larsen 
2009; see also Bland-Hawthorn et al. 2010b for a detailed dis- 
cussion). This mass range and value of 7 define the fiducial 
ICMF in the present investigation. 

The vast majority of star clusters in the Milky Way is found 
to be internally chemically homogeneous to a high degree (De 
Silva et al. 2006, 2007). This is to be expected from time 
scale arguments (Bland-Hawthorn et al. 2010b) and is incor- 
porated in the modeling by forcing each molecular cloud to 
be chemically homogeneous during its entire star formation 
period. This implies that a volume U c id = M/p, associated 
with each cluster, is made chemically homogeneous before 
star formation commences. Here, M = M*/e is the mass of 
the molecular cloud and p is the mean density of the inter- 
stellar medium (ISM). The star formation efficiency is set to 
e = 0.1. 

Associated with each cluster is also a mixing volume, V m i x , 
which is defined as the volume enriched by the ejecta from 
the SNe in the cluster. The ejecta are assumed to mix homo- 
geneously with the gas within V m i x . The size of the mixing 
volume can be expressed as (cf. Karlsson 2005; Karlsson et 
al. 2008) 

**( /3£^M r y /3 y /2 

1 / roix= _U ArbXt+ ^ j 1 , (2) 

where D tl h (kpc 2 Myr -1 ) denotes the turbulent diffusion co- 
efficient of the interstellar gas and M r denotes the typical 
mass swept up by a single SN remnant when the gas has 
cooled enough to form new stars. This occurs roughly when 
the remnant merges with the ambient medium, which implies 
that M r — 10 5 M Q (Ryan et al. 1996). We have introduced a 
small amount of randomness in the calculation of V m i X inso- 
far as log (M r ) is drawn from a normal distribution centered 
on log(Af r ) = 5 with a width of 0.25. For simplicity, we set 

Arb = 0. 

The formation of clusters is assumed to occur randomly in 
the simulation box. Each cluster, randomly distributed ac- 
cording to Eq. (Q~|i, is formed from a gas cloud with volume 
Vrfd which is made chemically homogeneous. If the cluster 
is massive enough to form high-mass stars (with masses dis- 
tributed according to the Salpeter IMF), the combined ejecta 
from the SNe are then mixed homogeneously within the vol- 
ume Vmix, given by Eq. (0. With the formation of each clus- 
ter, an amount of / out x M* of gas is subtracted from the 
galaxy's remaining gas mass. This includes mass expelled 
from the galaxy and mass locked up in long-lived stars and 
stellar remnants. The coefficient / out is empirically deter- 
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mined such that the peak of the model MDF is matched to 
the peak of the observed MDF (see Fig. |3J. For Sextans, 
/out = 37.4. And so, the chemical enrichment proceeds. Star 
formation is finally terminated when the gas density reaches 
p = 1 x 10 5 Mq, unless there is too little gas in the system 
for the final cluster to form. Further details of the model can 
be found in Bland-Hawthorn et al. (2010a). 

3.2. Stellar mass of Sextans 

We can estimate the present stellar mass of the Sextans dSph 
simply by adding up the light from its stars. The luminosity of 
Sextans is 4.1 ± 1.9 x 10 5 Cq, which corresponds to an abso- 
lute magnitude of My = — 9.2 ± 0.5 (Irwin & Hatzidimitriou 
1995). Since the majority of the stars in Sextans is believed 
to be older than 10 Gyr (e.g. Lee et al. 2009), we can esti- 
mate the mass simply by integrating the light from a single 
stellar population (SSP) of 12 Gyr, assuming a metallicity of 
Z = 10 _2 Z©. In this way, we approximate the present stellar 
mass to M Sx t = 8.9±4.1 x 10 5 Mq for a Salpeter IMF. Here, 
the symmetric uncertainty in luminosity is directly translated 
to an uncertainty in mass. 

We can, however, do slightly better and account for the ob- 
served spread in stellar age. Lee et al. (2009) determined the 
star formation history (SFH) of Sextans and noticed a spread 
in the data of at least 6 — 8 Gyr (excluding blue stragglers), 
even if most stars were recovered in the oldest age bins. From 
their inferred star formation rate (SFR; their case B), we adopt 
an SFH with four age bins; 15 - 13 Gyr, 13 - 11 Gyr, 11 - 9 
Gyr, and 9 - 7 Gyr, containing 66%, 25%, 7%, and 2% of 
the stars, respectively. The stars in each bin are assumed to be 
adequately described by an SSP of 14 Gyr, 12 Gyr, 10 Gyr, 
and 8 Gyr, respectively. The metallicity of the SSPs are all 
set to 1Q~ 2 Zq, except for the 8 Gyr SSP, for which a metal- 
licity of 1Q~ 15 Zq is assumed. For the light integration, the 
isochrones by Marigo et al. (2008) are adoptedQ From this 
modeling of the stellar content, we determine the present stel- 
lar mass of Sextans to Ms xt = 9.4 ± 4.4 x 10 5 Mq (Salpeter 
IMF). This is consistent with the mass Afsxt = 8.5 x 10 5 M ©, 
obtained by Woo et al. (2008). The stated uncertainty for- 
mally includes, apart from the observed error in the lumi- 
nosity, an estimate of the uncertainty from the inferred SFR 
(fSxt.SFR = 0.05 x 10 5 Mq) as derived from the data of 
Lee et al. (2009), as well as the intrinsic uncertainty due to 
the small number of luminous red giants in galaxies of this 
size (osxt.intr = 0.4 x 10 5 Mq). Evidently, these latter two 
uncertainties have a negligible impact on the total uncertainty. 

We note that any SFH inferred from photometric data may 
be marred with model dependent, systematic errors, e.g. from 
theoretical isochrone fitting, that are difficult to get a handle 
on and which add to the random errors in the SFR. The un- 
certainty crsxt.SFR could therefore be larger than the adopted 
value given above. From a set of completely random SFHs, 
the uncertainty in the mass determination is estimated to 
no more than csxt.SFR = 0.3 x 10 5 Mq, assuming that 
the youngest stars in Sextans are 8 Gyr. In the case of a 
Kroupa IMF (Kroupa 2001), we get a stellar mass today of 
5.4 ± 2.5 x 10 5 Mq. Note that the total mass (within 2.3 
kpc) of the dark-matter dominated Sextans is on the order of 
2 - 4 x 10 8 M Q (Battaglia et al. 201 1). 

In order to model the chemical evolution and reconstruct 
the galaxy with the correct number of clusters, we need to 

3 See http://stev.oapd.inaf.it/cgi-bin/cmd_2. 2 



estimate the initial mass of the stellar populations in Sex- 
tans. In a 12 Gyr population, for instance, only stars below 
m ~ 0.83 Mq are still present. Roughly 40-60% of the stel- 
lar mass is therefore gone (excluding stellar remnants, which 
have a negligible contribution to the total luminosity), de- 
pending on the IMF. This implies that, for a Salpeter (Kroupa) 
IMF, the birth mass of the stellar component of Sextans was 
Msxunit = 1-6 ± 0.8 (1.3 ± 0.6) x 10 6 Mq. This mass will 
be used below to estimate the age-weighted, initial mass of 
the st ar cl uster. 

In ij4.ll we will, along with the estimate of the most prob- 
able mass of the cluster, attempt to get a handle on the var- 
ious uncertainties involved. To this end, a detailed distribu- 
tion describing the uncertainty in the mass estimate of Sex- 
tans is desired. The luminosity of Sextans determined by Ir- 
win & Hatzidimitriou (1995) has a large uncertainty due to 
confusion between the dwarf and a background galaxy clus- 
ter. Consequently, a gaussian uncertainty distribution, implied 
by the stated symmetric error in luminosity, will have a non- 
negligible tail of negative luminosities. The probability of 
obtaining such false detections can, however, be suppressed 
by considering integral fluxes only over smaller regions cen- 
tered on the galactic core. As these fluxes are less affected 
by background noise, the probability of a non-detection is re- 
duced, keeping in mind that the mass inferred from these cen- 
tral fluxes is lower than the total stellar mass of the galaxy. 

Taking this suppression into account and retaining the stan- 
dard deviation in the luminosity stated by Irwin & Hatzidim- 
itriou (1995), we assume that the mass probability density 
function of Sextans is better approximated by a generalized 
Poisson distribution rather than the normal distribution. Its 
functional form is e~ mt (mt) kR /r(fcR + 1), where m de- 
notes the mass, t is a scale factor, and the parameter kji is 
a non-negative real number. Factorial k R is governed by the 
gamma function T(kji + 1). The specific values of and t 
are determined from the two relations k-^/t = 1.6 x 10 6 Mq 
and y/(k R + l)/t 2 = 0.8 x 10 6 Mq (Salpeter IMF), 
describing the most probable mass (mode) and the stan- 
dard deviation, respectively. It follows that /cr = 5.5 and 
t = 3.3 x 10~ 6 Mq" . As compared to a gaussian, this dis- 
tribution has a more extended tail towards high masses, dis- 
tinct from the suppressed tail towards low masses. It retains a 
fairly symmetric shape, nonetheless. Only a detailed study of 
how the background flux varies over the field in the immedi- 
ate vicinity of Sextans, can give us further information on the 
actual shape of the mass probability density function. 

4. NATURE OF CLUSTER 

We shall now turn to the relic star cluster itself. In the fol- 
lowing sections, we will estimate the initial stellar mass of 
the cluster and briefly discuss what possible type of cluster it 
might have been. 

4.1. Mass of cluster 

The mass of the dissolved cluster can be estimated indepen- 
dently from the size of the stellar "bump" in the MDF inferred 
from the re-calibrated DART data (Starkenburg et al. 2010) 
and from the identification of cluster stars in the data by Aoki 
et al. (2009). Let us first estimate the mass from the "bump" 
by comparing the simulated MDF with the observed one and 
then calculate the combined estimate, including the informa- 
tion on the num ber of "blue" stars identified as potential clus- 
ter stars in il2.ll 
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4.1.1. Mass estimate from the re-calibrated DART MDF 

As sugg ested by the appearance of the "blue" stars identified 
in jj2JJ we will assume that the "bump" predominantly origi- 
nates from a single, massive star cluster. This will have certain 
implications for the interpretation of the shape of the MDF. In 
the DART program, the [Fe/H] abundances are inferred from 
the CaT equivalent widths, which means that they de facto 
are translated [Ca/H] abundances. From Starkenburg et al. 
(2010, their Appendix A), it is clear that this translation sim- 
ply reads as an offset very close to 0.25 dex (i.e., [Ca/Fe] = 
+0.25±° l dex) over the entire range -4 < [Fe/H] < -0.5, 
for which the CaT calibration is valid. Any observed devi- 
ation from this presumed Ca-to-Fe enhancement (recall that 
the potential cluster stars have [Ca/Fe] ~ +0.05) would lead 
to an erroneous [Fe/H] determination. Hence, a translation 
like the one above merely shifts the DART MDF to the left in 
[Ca/H] -space, it does not reshape it into the form it actually 
should have if displayed in true [Fe/H] -space. In the mass- 
estimation analysis that follows below, we therefore choose 
to replace [Fe/H] with [Ca/H] as a measure of metallicity, in 
order to get a consistent comparison between models and ob- 
servations (see Fig. By applying a constant offset (suf- 
ficient for our purposes), the (re-)conversion to observed Ca 
abundances reads [Ca/H] b s = [Fc/H]dart + 0.25. 

The simulated MDF is composed of two main compo- 
nents; an underlying MDF of the galaxy and a cluster placed 
at [Fe/H] = -2.7 (i.e., [Ca/H] = -2.65, assuming that 
[Ca/Fe] = 0.05). A significant fraction of the stars in the 
sample of Aoki et al. (2009) are not identified with the rem- 
nant cluster, as evident from Fig. Q] The star least likely to 
be a cluster member (color coded red in Figs. Q]and[2]i is the 
most metal-poor star in the sample with [Fe/H] = —3.10. 
Moreover, the weighted mean metallicity of the "non-blue" 
stars is [Fe/H] = -2.90 ± 0.12, while that of the "blue" stars 
is [Fe/H] = -2.70 ± 0.12. This implies that the "bump" 
of cluster stars must lay on top of a broader distribution that 
extends down towards the extremely metal-poor regime (see 
Fig. [3]). We modeled the under lying MDF using the chemical 
evolution model described in ^3 - 1 1 assuming a Salpeter IMF 
and primordial (Z = 0), metal-independent SN yields from 
Nomoto et al. (2006). The Ca yields were divided by a factor 
of 1.8 {— —0.25 dex) to match the mean [Ca/Fe] observed in 
the Galactic halo. Similarly, the Mg yields were divided by a 
factor of 1.5 (= —0.18 dex). Since we are primarily interested 
in the chemical evolution below [Fe/H] ~ —2, the enrichment 
due to intermediate-mass stars and SNe type la were not in- 
cluded in the modeling. A total mass of ~ 1.6 x 10 6 Mq was 
imposed for the system "cluster + underlying galaxy". 

By matching the number of observed stars in the "cluster 
region", arbitrarily defined as the region in the metallicity in- 
terval -3.15 < [Ca/H] < -2.15, with the model MDF (inte- 
grated over the same range), a fiducial cluster mass was deter- 
mined. The procedure is illustrated in Fig. [3] The integrated 
number of observed stars within this region is more than 3a 
above the model of the underlying galaxy while the number of 
stars in the central metallicity bin at [Ca/H] = —2.65 is more 
than 4er above the underlying galaxy. The "bump" (light-gray 
shaded area) contains 12.3% of the total stellar mass, which 
implies an initial cluster mass of M* i i n it = 2.0 x 10 5 Mq. 
For a Kroupa IMF, the initial cluster mass is slightly lower, 
or M^init = 1.6 x 10 5 Mq. These masses translate into 
present-day cluster masses of M* = 1.2 (0.7) x 10 5 Mq, for 
a Salpeter (Kroupa) IMF. 




Figure 3. Metallicity distribution function for stars in the Sextans dSph 
galaxy, as given by d/»/d[Ca/H]. The quantity /„ is the fraction of stars 
that fall below each [Ca/H] bin (1 dex). The black, solid line denotes the fidu- 
cial distribution of [Ca/H], predicted from the stochastic chemical evolution 
model presented in ij3]while the black, dashed line denotes the corresponding 
distribution, including an Af* = 2.0 X 10 Mq star cluster at an assumed 
Ca abundance of [Ca/H] = —2.65 (arrow). Both MDFs are convolved with 
a a = 0.2 dex gaussian. The red step function (with Poissonian noise) shows 
the observed distribution of [Ca/H], as inferred from the DART data (Starken- 
burg et al. 2010). Note that the data point at [Ca/H] = -2.65 lies more than 
4<r above the model of the underlying galaxy. 

In order to get a proper handle on the uncertainties involved 
and to estimate the most probable cluster mass, we must ac- 
count simultaneously for the uncertainties of all parameters 
that go into the calculation of the cluster mass. Given the 
observed and modeled MDF, the mass of the star cluster is 
determined from the expression: 



M, 



Mi 



Sxt,init 



4 + iv Q 



(3) 



where 7Vj n denotes the observed number of stars within the 
"cluster region", N out denotes the observed number of stars 
outside this interval, and Cent denotes the continuum level 
defined as the mean number of stars (normalized to the ob- 
servations) in the part of the underlying model galaxy that 
falls within the "cluster region". The mass of Sextans is de- 
noted by A/sxt.init- Note that all parameters in Eq. (f3]l are 
random variables (here, explicit random variables are indi- 
cated by the "hat" symbol), each to which a specific distri- 
bution must be specified in order to determine the distribution 

of cluster masses, M*,i n it. 
To evaluate the uncertainty in the cluster mass associated 

with the observations of the MDF, the random variables N- ln 
and iVout are assumed to be distributed according to the Pois- 
son distribution e ^fi^/kl, with /i; n = 71 and /i ou t = 111, 
respectively. These numbers are extracted from the DART 
data of Sextans (Starkenburg et al. 2010). The correspond- 
ing dispersion in the cluster mass due to this uncertainty is 
c*,init = 0.7 x 10 5 Mq. Moreover, as a result of stars being 
formed in clusters, the shape of the model MDF of the un- 
derlying galaxy should vary somewhat from galaxy to galaxy. 
We collected 100 di ffere nt realizations, assuming the fiducial 
ICMF discussed in £13.11 to map out this effect. 

With the assumption that the "bump" at [Fe/H] = —2.7 
is predominantly due to a single cluster, we accepted only 
relatively smooth shapes, disregarding any model with obvi- 
ous "secondary bumps" in the metallicity range of interest. 
The variations in the continuum level were found to be nearly 
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Poisson distributed around the mode (i.e., the continuum level 
of the fiducial model depicted in Fig. with a standard de- 
viation of Cent = 6.2. We note that with such a a, the dis- 
tribution is fairly close to being normal. It follows that the 
uncertainty of the cluster mass due to the underlying model 
galaxy alone is o"* fait = 0.7 x 10 5 Mq. This uncertainty 
depends on parameters like the ICMF which have an impact 
on the shape of the MDF. A larger fraction of smaller clusters 
would, for example, result in a smaller dispersion of the clus- 
ter mass. Note that, for an ICMF with a slope departing from 
7 = 2, the fluctuations in the continuum level may be de- 
scribed by a distribution that deviates from that of the Poisson 
distribution. A strongly varying star formation rate would also 
affect the continuum level, but this variation is not included in 
the present investigation. Finally, from the uncertainty in the 
mass of Sextans, governed by the distribution introduced in 
A3. 21 the uncertainty in the mass of the cluster is estimated to 
ff*,teit = 0.9 x 10 5 Mq. 

When accounting for the uncertainties in the observed and 
modelled MDF and the mass of Sextans simultaneously, we 
derive, from the expression in Eq. ([3), a most probable clus- 
ter mass of M*,fa it = 1.6±^ x 10 5 Mq (Salpeter IMF). This 
mass is slightly lower than the fiducial mass deduced above. 
The asymmetric uncertainties are calculated such that there is 
a 68.3% probability of having a cluster mass within the stated 
interval, where the endpoints denote the intersections between 
the distribution and a horizontal line (cf. Fig. @). Note that 
the black, dashed line in Fig. |4]exhibits a non-zero probabil- 
ity having negative cluster masses. This just reflects the small 
statistical chance that the continuum level of the underlying 
galaxy may lie above the observed MDF. In these cases, no 
cluster is detected. If we take this into account, the mean clus- 
ter mass is 2.3 x 10 5 Mq while the median is 2.1 x 10 5 Mq. 

4.1.2. Combined mass estimate 

We note that the predicted fraction of cluster stars below 
[Fe/H] = —2.5 in the particular model displayed in Fig. [3] 
is 40%. This implies a nearly l-in-2 chance (binomial statis- 
tics) of finding three, or more, cluster stars in a random sam- 
ple of six stars in this metallicity range. In less than 5% of 
the cases, zero cluster stars are expected in a sample of six 
stars. Clearly, this is consistent with the recognition of th ree 
(or four, including S 11 — 37) potential cluster stars ( fl2.ll ) in 
the data of Aoki et al. (2009). The likelihood of this indepen- 
dent finding may be utilized to put additional constraints on 
the mode and mass range of the cluster, as shown in Fig. [4] 
It should be pointed out that the stars in the sample of Aoki 
et al. (2009) only probes a fraction of the entire stellar pop- 
ulation of Sextans and the local fraction of "blue" stars may 
very well change with distance from the galactic center. As 
a zeroth-order approximation, we will, however, assume that 
the fraction of former member stars of the relic cluster at a 
distance 0°.4 < R < 0°.6 is equal to the average fraction 
taken over the entire galaxy out to R = 1°.8, which is the 
outer boundary in the DART survey (Battaglia et al. 201 1). 

At face value, the identification of the three "blue" stars in 
Fig. Q] if they do belong to a dispersed cluster, point towards 
a cluster mass on the order of a few hundred thousand so- 
lar masses. The probability that exactly three out of six stars 
below [Fe/H] = —2.5 should belong to the cluster is higher 
than for any other choice of k (i.e., assuming a Salpeter IMF 
and a galaxy mass of Mgxt.fait = 1-6 X 10 6 Mq) in the mass 
range 2.2 x 10 5 < M*,fa it /!M < 3.6 x 10 5 . That is to say, in 
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Cluster mass, M^ irlit (x 10 5 Mq) 

Figure 4. The uncertainty in the estimation of the cluster mass. The black, 
dashed line denotes the distribution derived from Eq. while the black, 
solid line denotes the combined, conditional distribution given that three out 
of six stars below [Fe/H] = —2.5 are identified as cluster stars in the data of 
Aoki et al. (2009). The probability for this to occur is indicated by the gray, 
dashed line. The most probable mass is determined to 1.9 g x 10 5 Mq for 
a Salpeter IMF. The (asymmetric) uncertainty in the estimation is indicated 
by the horizontal arrows. There is a 68.3% probability that the cluster mass 
falls within this interval (shaded area). 

this mass range, the binomial distribution B'm(pM, ,k,n = 6) 
peaks at k = 3, where p M , = M*,init/Msxt,init | [Fe/H]<-2.5- 
Taking into account the uncertainty in the mass of the galaxy, 
the most probable cluster mass is M*,jnit = 3.1 x 10 5 Mq for 
k = 3. Given this conditional constraint, the combined esti- 
mate of the most probable, initial mass of the dissolved clus- 
ter (with asymmtric uncertainties, see Fig. |4]( is determined 
to M^fait = 1.9±o;l (1-6±J:|) x 10 5 Mq, for a Salpeter 
(Kroupa) IMF. This is our most reliable estimate of the mass 
of the relic star cluster. 

4.2. Type of cluster 

Besides its relatively large initial mass, can we say some- 
thing more specific about the possible nature of the relic star 
cluster? Most Galactic open clusters are found in the mass 
range 10 2 < M^/Mq < 10 4 , but more massive, young star 
clusters do exist, like Westerlund 1. This cluster has a dy- 
namical mass of 1.5 x 10 5 Mq (Mengel & Tacconi-Garman 
2009). Another, intermediate-age cluster, namely GLIMPSE- 
CO 1, has an estimated mass of 8 x 1O 4 A^ (Davies et al. 
2011). These clusters are roughly in the same mass range as 
the remnant cluster in Sextans. Moreover, even if the aver- 
age present-day mass of the open clusters in the Milky Way is 
about 7 x 10 2 Mq, the average mass of the ICMF could have 
been as high as 4.5 x 10 3 Mq, owing to efficient evaporation 
(Piskunov et al. 2008). Although rare, open clusters with ini- 
tial masses on the order of 10 5 Mq are therefore conceivable, 
at least in the Galactic Disk. 

The masses of Galactic globular clusters (GCs, hereafter) 
are, on the other hand, found in the range 1 x 10 4 < 
M*/M Q < 3 x 10 6 . The Sextans cluster falls right in the 
middle of this range. GCs also show specific elemental abun- 
dance correlations, such as the Na-O and the Al-Mg anticor- 
relations. It has recently been established that these anticorre- 
lations are caused by internal pollution of the intracluster gas 
(see Carretta et al. 2009 for a discussion). Such correlations 
are not observed in open clusters and could therefore, hypo- 
thetically be used to discern the nature of the relic cluster. Un- 
fortunately, Aoki et al. (2009) were neither able to measure 
oxygen, nor aluminum in any of the stars in their sample. 

Interestingly, however, the "blue" star S 10—14 has a signif- 
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icantly lower Na abundance than the other two potential clus- 
ter stars, with a measured Na-to-Fe ratio of [Na/Fe] = —0.72 
(see Fig. [TJ. In comparison, S 11 — 13 and S 14 — 98 both 
have [Na/Fe] ~ —0.1. This may be suggestive of an Na-O 
anticorrelation present in the data. We note that all the "blue" 
stars are found at the lower end of the [Na/Fe] -distribution 
(i.e., LTE abundances) usually observed in globular clusters. 
Aoki et al. (2009) interpreted the absence of Na enhance- 
ments as evidence against a GC origin for the Mg deficient 
stars. However, if the primordial Na abundance was low in 
the gas cloud out of which the GC was formed, low Na abun- 
dances on average could be expected. Such an effect is also 
present in the sample of Galactic GCs observed by Carretta et 
al. (2009). Furthermore, Carretta et al. (2009) conclude that 
roughly one third of the stars in Galactic GCs exhibit primor- 
dial abundances, while two thirds are affected by internal pol- 
lution by the GCs' first stellar generation. This is consistent 
with the data of Aoki et al. (2009), assuming that S 10 - 14 
belongs to the first, primordial generation while S 11 — 13 and 
S 14 — 98 belong to the second, polluted generation. 

While the vast majority of the Galactic GCs shows no ob- 
served variations in the neutron capture elements between in- 
dividual member stars, two of the most metal -poor GCs, Ml 5 
(Sobeck et al. 2011) and M92 (Roederer & Sneden 2011), 
have recently been demonstrated to exhibit a spread in a num- 
ber of heavy elements, including Ba and Eu. The mechanism 
behind this star-to-star scatter is not known. We note that at 
least 80% (excluding upper limits) of the stars in the samples 
of Sobeck et al. (2011) and Roederer & Sneden (2011) are 
r-process enhanced (r-I) and all stars have [Ba/Fe] > —0.25. 
In contrast, the mean ratio of the potential cluster stars is 
([Ba/Fe]) = -1.33. Moreover, both M15 and M92 are sev- 
eral times more massive than what the relic cluster would have 
been, if it had survived to the present epoch. 

It is unclear whether all GCs below [Fe/H] ~ —2.3 should 
exhibit chemical inhomogeneities in the neutron capture ele- 
ments. Letarte et al. (2006) observed three very metal-poor 
GCs in the Fornax dSph, of which two were found to have 
[Fe/H] < —2.3. Their data show no evidence in support of 
an intrinsic scatter in any of the GCs, although the statistics 
is admittedly low. At this point, we can therefore neither rule 
out the possibility that the relic cluster was a GC based on 
the observed clumping of stars in [Ba/Fe], nor that additional 
stars, e.g. S 12 — 28 (color coded yellow in Fig. [TJ, should 
be included as potential members. Further observations, in- 
cluding measurements of O, Al, and neutron capture element 
abundances are required in order to resolve these issues and 
determine the true nature of the relic cluster. 

4.3. Possible origins of the atypical [a/Fe] ratio? 

In Fig. [TJ the [Mg/Fe] and [Ca/Fe] ratios of the "blue" 
stars, as well as of S 11 — 37, are unusually low for stars 
in this metallicity regime (e.g. Tolstoy et al. 2009, their 
Fig. 16). In the present modeling, we simply appended an 
-M*,init = 1-9 x 10 5 A4q star clust er at t he observed loca- 
tion of the "blue" stars in C-space ( 34.1.11 see also Fig. 0. 
However, the probability for a cluster of this mass to be found 
in this part of C-space (i.e., in [Ca/Fe]) purely by chance is 
estimated to < 1%, given the fiducial model and that nor- 
mal core-collapse SNe were the only nucleosynthesis sites 
enriching the early ISM. We note in passing that relatively 
low [Ca/Fe] ratios are found in five Galactic Halo GCs (see 
compilation by Pritzl et al. 2005). Of these, only the two 
most metal-poor clusters; NGC 2419 ([Fe/H] = -2.32) and 



NGC 5466 ([Fe/H] = -2.05, both have [Ca/Fe] = +0.1) 
are believed not to have been captured from a satellite system 
with a slower chemical evolution. The Mg-to-Fe ratio of the 
two clusters is determined to be +0.30 and +0.06, respec- 
tively. Unfortunately, the abundance analysis in both cases is 
based on one star only. 

A decrease of the [a/Fe] ratio with increasing metallicity is 
commonly believed to be a result of the delayed enrichment 
by SNe la synthesizing iron and Fe-peak elements (Tinsley 
1980; Matteucci & Greggio 1986). Due to the slower star for- 
mation timescale in dSphs, a subsolar [a/Fe] is reached at a 
lower [Fe/H], as compared to the Galactic Disk. However, in 
general, the classical dSphs only exhibit stars with a subsolar 
[a/Fe] ratio above [Fe/H] ~ -1.5 (Tolstoy et al. 2009; Kirby 
et al. 2011a). At [Fe/H] ~ —3, typical [a/Fe] enhancements 
are observed, also in the dSphs (see Fig. [2J. In most ultra- 
faint dwarf galaxies, we have a similar situation with the bulk 
of stars being [a/Fe] -enhanced below [Fe/H] ~ — 2 (Feltzing 
et al. 2009; Frebel et al. 2010b; Norris et al. 2010a, c; see, 
however, Aden et al. 201 1). 

In a dwarf galaxy with a very low SFR, it is conceiv- 
able that SNe la start contributing to the metal budget below 
[Fe/H] = -2 (see Revaz et al. 2009). Alternatively, infall of 
metal-free gas could in principle lower the [Fe/H] of an SN 
la-dominated ISM to the desired metallicity regime. Ivans et 
al. (2003) reported on three Galactic halo stars at a metallicity 
[Fe/H] ~ — 2 with [a/Fe] ratios similar to those observed in 
the "blue" stars. They concluded that the abundance patterns 
of these stars are consistent with an early enrichment by SNe 
la. But on comparable timescales as the SNe la, the asymp- 
totic giant branch (AGB) stars would also contribute to the 
metal budget. These stars enrich the ISM mostly in C, N, and 
s-process elements like Ba. Aoki et al. (2009) therefore re- 
ject the SNe la as being responsible for the low [a/Fe] ratios 
due to the low [Ba/Fe] ratio observed in the "blue" stars. Fur- 
thermore, in a medium-resolution spectroscopic study, Kirby 
et al. (201 la) mapped out the chemical evolution of 8 dSphs, 
including Sextans. In their data for Sextans, [Mg/Fe] reaches 
subsolar values at [Fe/H] ~ —1.6, while [Ca/Fe] ~ is 
reached at somewhat lower [Fe/H]. They point out, however, 
that the [a/Fe] ratio appears to be affected by SNe la already 
at [Fe/H] ~ —2.5. At such a low metallicity, small num- 
ber statistics plays a considerably more prominent role than 
it does at [Fe/H] > —1.5. The parent molecular cloud from 
which the cluster was born could theoretically have been en- 
riched in SN la material, but not much in AGB material, as 
implied by Ivans et al. (2003) for their Galactic halo stars. 
Such a scenario must be checked against detailed modeling of 
the early chemical evolution of Sextans. We plan to do so in a 
forthcoming paper. 

As suggested by Aoki et al. (2009), a truncated IMF at the 
highest masses would favor lower [a/Fe] ratios in the star- 
forming gas. Weidner & Kroupa (2005; see also Pflamm- 
Altenburg et al. 2007) argue that systems with slow star for- 
mation will effectively acquire a truncated IMF as a result 
of the relation between the SFR and the maximum cluster 
mass, given that the mass of the most massive star formed 
in a cluster is limited by the mass of the cluster. Since [a/Fe]- 
enhanced stars do exist in the sample (S 15 — 19), some high- 
mass SNe must, nonetheless, have contributed to the metal 
budget in Sextans. A similar result is obtained for an IMF 
with a steeper than usual slope at the high-mass end. 

Could the presence of hypernovae (HNe) contribute to a 
larger fraction of solar- [a/Fe] stars? As higher temperatures 
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Figure 5. Predicted distributions of stars in C-space for a Sextans-sized model galaxy. The top panels show the probability of finding a star in the [Mg/Fe]-[Fe/H] 
plane (left) and the [Ca/Fe]-[Fe/H] plane (right). The big, blue circle indicates the location of former member stars of a M» = 1.9 X 10 5 Mq relic cluster, while 
the small, red circles indicate the locations of stars which once belonged to clusters that were formed out of gas only enriched by a single, primordial SN, given 
that the first stars in dSphs were formed from primordial gas. The predicted fraction of these second generation stars in galaxies like Sextans is 1.1 X 10~ 3 , while 
the fraction below [Fc/H] = —2.5 is 3.7 X 10 , for the fiducial ICMF. Due to the stochastic nature of the chemical enrichment, second generation stars do not 
necessarily show up in the most metal-poor regime at [Fe/H] < —4. An intrinsic uncertainty of no more than 0.01 dex is applied to all abundance ratios in order 
to emphasize the presence of clumps. The bottom panels show the same distributions convolved with a two-dimensional gaussian with a = 0.15, 0.20, 0.20 
dex in [Mg/Fe], [Ca/Fe], and [Fe/H], respectively. Most of the clumping is washed out, except the one originating from the "blue" star cluster. The size (la- 
uncertainty) and location of this clump is indicated by the dash-dotted, blue circles. White circles denote the six stars observed by Aoki et al. (2009), while the 
squares and diamonds denote observations by Tafelmeyer et al. (2010), and Shetrone et al. (2001), respectively. The six big, colored circles are model stars 
randomly drawn from the non-shaded part (i.e., [Fc/H] < —2.5) of the distribution functions in the top panels. The three stars displayed in blue originate from 
the M* = 1.9 X 10 Mq cluster, while the black stars originate from other clusters. Fiducial uncertainties of 0.15 dex in [Mg/Fe], 0.20 dex in [Ca/Fe], and 
0.20 dex in [Fe/H] have been applied, similar to those reported by Aoki et al. (2009). 



are realized in models of HNe (Umeda & Nomoto 2002; 
Nomoto et al. 2006), the complete Si-burning region ex- 
tends outwards in mass, which leads to larger amounts of Fe 
( 56 Ni) in the ejecta, while the light Fe-peak elements Cr and 
Mn are underproduced. Lower [a/Fe] ratios can therefore be 
achieved for the high-mass SNe, even though the lowest ratios 
are still maintained in the low-mass SNe. The larger explo- 
sion energies associated with the HNe also result in the ejecta 
being diluted with a larger mass of ambient gas before merg- 
ing with the ISM. On a local scale, the presence of HNe thus 
induce more efficient mixing which decreases the possibility 
of having star-forming regions with [a/Fe] ~ at any given 
[Fe/H]. This implies that even fewer solar- [a/Fe] stars should 
be able to form around [Fc/H] = —2.7, despite lower [a/Fe] 
ratios in the ejecta of these highly explosive SNe. Alterna- 
tively, if the HN ejecta are lost from the system as a result of 
the higher exp losions energies (along with a shallow potential 
well, see ^2. 3b . a lower [a/Fe] ratio would also be realized in 
the remaining, mixed gas, although the overall effect would 
be smaller. Consistent with the HN scenario, a relatively low 
[Cr/Fe] ratio is measured in the "blue" stars (Fig. [2]), as com- 
pared to the bulk of dSph stars. Why HN enrichment then 
appears not to have occurred in all dSphs should be further 
investigated but could be explained by local, chemical inho- 
mogeneities in the star-forming gas. 

There is an interesting alternative explanation to the ob- 
served [Mg/Fe] and [Ca/Fe] ratios in the "blue" stars. These 



stars could have formed out of gas that was partly enriched by 
a very massive Pop III star exploding as a pair-instability SN 
(Heger & Woosley 2002; Umeda & Nomoto 2002). Yields 
for the most massive stars in the pair-instability mass range 
(i.e., 240 < m/M.® < 260) produce roughly solar and 
subsolar [Ca/Fe] and [Mg/Fe] ratios, respectively (Heger & 
Woosley 2002). Karlsson et al. (2008) predict that stars with 
a predominant contribution from primordial pair-instability 
SNe should have relatively high Ca abundances roughly in 
the range —3 < [Ca/H] < —2, in contrast to what may 
be expected for stars enriched by a first stellar generation. 
This result is also supported by hydrodynamics simulations 
which obtain an overall metallicity of Z ~ 10~ 3 Z in proto- 
galaxies enriched by pair-instability SNe (Wise & Abel 2008; 
Greif et al. 2010). Evidently, the Ca abundance of the "blue" 
stars is found within the predicted range. Furthermore, as a 
result of a very small neutron excess in the interior layers, 
low [Na/Fe], [Sc/Fe], and [Ba/Fe] ratios (see Fig. [TJ are to be 
expected the ejecta of pair-instability SNe. The pronounced 
odd-even effect is particularly interesting for Na since it could 
explain the very low Sodium abundance in S 10 — 14. 

Further observations including measurements of C, O, Si, 
and Zn abundances are of crucial importance in order to dis- 
criminate between the various scenarios presented here. Ad- 
ditional information on stellar ages would also help us to de- 
termine whether gas infall has undermined a closed-box-like 
chemical evolution of Sextans (see, e.g. the discussion on 
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SNe la). As a result of the very high explosion energies in 
pair-instability SNe (~ 10 52 — 10 53 erg), the oxygen-burning 
region in the models is significantly extended as compared to 
the corresponding region in core-collapse SNe, and even HNe. 
Oxygen-burning products like Si and S should therefore be 
abundant in the ejecta relative to C and O and high ratios of, 
e.g. [Si/O] and [S/C] would be expected. As no fall-back onto 
a remnant occurs (the stellar core is complete disrupted in the 
explosion), the [O/Fe] ratio is predicted to be subsolar in stars 
predominantly enriched by pair-instability SNe while it is su- 
persolar in stars enriched by core-collapse SNe and HNe, in 
line with with the bulk of the stars in the metal-poor Galactic 
Halo. Zn may be another discriminator, which is predicted to 
be produced in abundance in HNe, but not in pair-instability 
SNe (Umeda & Nomoto 2002). 

4.4. Remark on the cluster signature in C-space 

Figure[5]illustrates the expected imprint of a 1.9 x 10 5 A4q 
cluster in C-space, assuming that the early star formation in 
Sextans was clustered. The upper panels indicate the rich 
spectrum of clumps that is predicted to exist. This is in con- 
trast to predictions in which clusters either are not accounted 
for (e.g. Karlsson & Gustafsson 2005; Cescutti 2008) or the 
mass distribution of the clusters is not accounted for (e.g. Re- 
vaz et al. 2009). These models produce smoother distribu- 
tions in abundance space, although the simulations by Revaz 
et al. (2009) exhibit some clumping. However, due to the 
observational uncertainties typical for present-day studies of 
extra-galactic stars, the wealth of structure in the upper pan- 
els of Fig. [5] will almost entirely be washed out, except for 
the 1.9 x 10 5 M.q cluster (lower panels). Uncertainties of 
< 0.1 dex must be achieved in order to reach the level of ac- 
curacy required for multiple clumping to appear in a-element 
abundance space (Bland-Hawhorn et al. 2010; Karlsson et al. 
201 1). The prediction that underlying clumping exists below 
the current detection limit is a key point. 

Note that the three blue model stars in the lower panels of 
Fig. [5] clump rather tightly together in Mg but slightly less 
so in Ca. A situation like this may occur purely by chance 
due to the observational uncertainties and illustrates the need 
for studies, such as the one by Aoki et al. (2009), in which 
accurate stellar abundances of many elements are measured. 
Having access to a multidimensional C-space in combination 
with a large stellar sample also minimizes the risk of detecting 
"false clumps" in individual elements, which otherwise could 
lead to erroneous conclusions about the stars' origins. 

But if we do not have access to accurately determined abun- 
dances of many elements, can the cluster signature be detected 
in one dimensional C-space? We shall claim that this is the 
case, at least in a statistical sense, and argue that the cluster 
signature may already be prevalent in the distributions of stel- 
lar metallicities of dSphs (see fJ5]below). 

5. DISCUSSION AND IMPLICATIONS 

The probable identification of stellar debris from a dissolved, 
ancient cluster at [Fe/H] = —2.7 in the Sextans dSph have 
several interesting implications for the understanding of early 
star and galaxy formation. In the following section, we shall 
discuss the interpretation of MDFs and cumulative metallic- 
ity functions of dwarf galaxies from a statistical point of view, 
particularly in connection to the formation of star clusters. We 
will also briefly discuss what implications this has for the un- 
derstanding of star formation in metal-poor systems, for the 



origins of dwarf galaxies, and for the build-up of the stellar 
Halo around the Milky Way. 

5.1. Implications for early star and cluster formation 

We know little about the star formation process at the earliest 
cosmic times. Inefficient fragmentation of collapsing, metal- 
free gas clouds was believed to prevent the formation of pri- 
mordial clusters (Yoshida et al. 2008). However, in more 
recent hydrodynamics simulations, primordial stars are found 
to be formed in binary or small multiple systems (Turk et al. 
2009; Stacy et al. 2010), although these simulations have not 
yet reached their asymptotic end state. Assuming turbulent 
initial conditions, Clark et al. (201 1) showed that primordial 
stars might even have been able to form in small, dense clus- 
ters. In their simulations, the resulting stellar IMF is also sig- 
nificantly broader than what previously have been found (e.g. 
Bromm et al. 2009; see, however, Nakamura & Umemura 
2001), with a noticeable fraction of stars below 1 A4q. See 
Karlsson et al. (201 1) for a recent review on the topic. 

Our knowledge of clustered star formation in the metallic - 
ity regime — 5 < [Fe/H] < — 2 is perhaps even more limited; 
here only a few of the most metal-poor Galactic globular clus- 
ters are found. Wise et al. (2012b) introduce Population II 
star clusters in their simulations of the birth of the first galax- 
ies, without actually discussing the cluster formation process 
itself. From their data (J. Wise, priv. comm.), we find that 
most clusters below [Z/H] = —2 were formed with masses 
in the range 7.8 x 10 2 < M*/M@ < 1.4 x 10 4 . The clus- 
ters below the low mass limit were not properly resolved and 
were rejected. The slope of the numerical ICMF in the above 
mass range is close to 7 = 2, and compatible with the slope 
of the present-day ICMF. This is, indeed, an interesting re- 
sult and predicts that 7 may be universal all the way down to 
[Fe/H] ~ —3.5, at least in certain environments. By analyz- 
ing the metal content of the stellar populations in old, metal- 
poor dwarf galaxies, we should be able to put constraints on 
these numerical results and push our knowledge of the star 
formation process into the most metal-poor regime. In the 
analysis that follows, we will, for convenience, adopt a mass 
range of 5 x 10 2 < M*/Mq < 1 x 10 4 and a slope of 7 = 2 
for this ICMF. 

If the i nterp retation of the chemical abundance data pre- 
sented is £12.1 [ indeed is correct, the "blue" stars in Fig. Q]were 
once members of the most metal-poor star cluster known to 
have formed in the Universe. This clearly suggests that stars 
could form in massive clusters already around [Fe/H] ~ —3. 
However, the question remains whether clustered star forma- 
tion was common already at these metallicities and whether 
the mass range and relative frequency of clusters resembled 
those of the present-day ICMF. 

Consider the MDF, or the cumulative metallicity function, 
/*. Suppose that star formation in the early epochs was highly 
unclustered, i.e., stars were predominantly forming in isola- 
tion. If so, the MDF, as well as /», would appear the same 
for all galaxies with identical (same SFR, infall/outflow rate 
etc.) evolution. However, the situation would be different if 
stars were formed in chemically homogeneous clusters. In 
this case, each MDF would exhibit a unique, irregular and 
"bumpy" shape as a result of stars grouping together in ran- 
dom clumps in metallicity-space. As an effect of small num- 
ber statistics, these irregularities should be most notable for 
the least massive galaxies, and at the lowest metallicities. 

We should emphasize that not all "bumps" must necessarily 
be signatures of disrupted clusters. Star bursts, gas infall, the 
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Figure 6. Predicted cumulative metallicity functions, /» (normalized to unity at [Fc/H] = —2.5), for Sextans-sized dSph galaxies. Each of the left-hand-side 
panels show 10 model /* (colored squares) for different values of the slope of the ICMF, for which 7 = 2 is the fiducial value. The top right panel shows the 
resulting cumulative functions assuming the narrower cluster mass range inferred from the simulations by Wise et al. (2012b), while the middle right panel shows 
the /* in the case of isolated, unclustered star formation. Each model /* is convolved with a a = 0.2 dex gaussian and sampled with 60 stars. The bottom right 
panel displays the observed /* (Starkenburg et al. 2010) for Carina (gray circles), Sculptor (squares), Sextans (down-faced triangles), and Fornax (left-faced 
triangles). The combined /* for the 8 ultra-faint dwarf galaxies (small, green diamonds) discussed by Kirby et al. (2008, 2010, 2011b) is plotted for reference. 
In all panels, the dashed line denotes the predicted mean, /», for a large number of model dSphs, while the black diamonds denote the corrected, cumulative 
function, /*,Halo (binned), observed for the Galactic halo (Schorck et al. 2009). The offset between /» and /*,Hak> is predominantly due to the lack of stars 
below [Fe/H] ~ —3.5 in the medium-resolution (1Z ~ 2, 000) HES data (see text). The sand-colored areas (bounded by dotted lines) are measures of the 
dispersion expected from the finite sampling (i.e., 60 stars) of the cumulative functions and denote the regions within which 80% of the /* are confined, given 
that they are drawn directly from /„ (light-shaded) or /»,Halo (medium-shaded, displayed for reference). The overlap between these two areas is shaded dark. 
In cases of unclustered star formation (middle right panel) or when the clusters are mostly "unresolved" (middle left and top right panels), the vast majority of 
the model /* is expected to be found within or in close connection to the light-colored region. 



onset of SNe la (see Karlsson 2005, their Fig. 10), and even 
statistical uncertainties in the observations could in principle 
introduce irregularities. However, "bumps" originating from 
variations in the global SFR or infall of gas may have notably 
larger intrinsic widths in metallicity space since they do not 
originate from a single, chemically homogeneous stellar pop- 
ulation. "SN la bumps" in [Fe/H] are formed when the first 
stars enriched by SNe la are pushed towards higher [Fe/H], 
relative to stars only enriched by core collapse SNe. This 
actually generates a through in the MDF, whose metal-poor 
wall appears as a "bump". Features like this will only show 



up in elements predominantly synthesized in the SNe la. In 
elements like Mg or Ca, the "SN la bump" is significantly 
diminished, in contrast to "bumps" originating from clusters. 
We will therefore argue that MDFs and the corresponding /* 
can be used to probe the early formation of clusters in a sta- 
tistical sense. 

5.1.1. Star formation in dwarf spheroidal galaxies 

We ran a series of simulations of Sextans-sized dSphs, assum- 
ing different ICMFs, to map out the effect. For simplicity, 
Sextans was chosen as a typical dSph in terms of luminosity 
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and stellar mass. In the sample of Starkenburg et al. (2010), 
Carina, Sextans, and Sculptor all have similar masses (within 
a factor of 2, see Woo et al. 2008). Only Fornax is signifi- 
cantly more massive, with a stellar mass roughly 20 x that of 
Sextans. These four dSphs will be used in the comparison to 
the simulations. Here, we are only interested in irregularities 
and "bumps" of the /„ (and the corresponding MDFs) in a sta- 
tistical sense. We will therefore take the observational data as 
given by Starkenburg et al. (2 010). N o translation to [Ca/H] 
abundances is performed (cf. 34.1 .11. The same goes for the 
ultra-faint dwarfgalaxy data in i )5.1.2l The results are sum- 
marized in Fig. [6] Each cumulative metallicity function, /», 
is convolved with a a = 0.2 dex gaussian and sampled with 
60 stars. The finite sampling is chosen to match the sampling 
of /» for Sextans (Starkenburg et al. 2010) and introduces 
an observational scatter indicated by the sand-colored areas 
(light-colored, see caption of Fig. [6] for details). The disper- 
sion seen in the middle right panel of Fig. [6] (unclustered star 
formation), is solely due to the finite sampling of the cumula- 
tive functions. 

At face value, the ICMF with a shallow slope of 7 = 1 (Fig. 
[6j bottom left) overproduces the number of model galaxies 
with a highly irregular /* as compared to observations (bot- 
tom right). In contrast, the ICMF with 7 = 2.5 (middle left) 
and the ICMF inferred from the simulation data of Wise et 
al. (2012b, top right) produce very few /* that significantly 
deviate from what is expected from the finite sampling alone, 
and the plots are virtually indistinguishable from that of the 
unclustered star formation (middle right). Higher sampling, 
i.e., more stars, and more precise metallicities are necessary 
in order to distinguish between these latter two ICMFs. 

The fiducial ICMF (top left) is an intermediate case. It leads 
to slightly pronounced fluctuations as compared to the case of 
unclustered star formation. In fact, about 75% of the /» in the 
top left panel is predicted to fall outside the boundaries (in this 
case defined by the light-shaded, sand-colored area) at some 
point in metallicity. For unclustered star formation (i.e., the 
scatter is solely due to statistical fluctuations) the correspond- 
ing fraction is 20%. Note that the sand-colored areas in Fig. [6] 
are defined such that 80% of all /*, drawn from a single parent 
distribution and sampled with a given number of stars (60 in 
the case of Sextans), should be confined within the limits over 
the entire metallicity range [Fc/H] < —2.5. The scatter gen- 
erated by the fiducial ICMF is in qualitative agreement with 
the observations of the four classical dSphs (Fig. [6] bottom 
right panel). At first glance, the scatter in the observations 
may look similar to the unclustered case. But a careful in- 
spection of the dSph cumulative metallicity functions reveals 
that all four /* at some point are found outside both the light- 
and dark-shaded, sand-colored areas, taking into account the 
poorer sampling of the /* of Sculptor, Carina and Fornax. As- 
suming unclustered star formation in the dSphs, the chance is 
less than 0.2% that four out of four (identical) /* should be 
found outside these limits. The corresponding probability as- 
suming a present-day ICMF in the dSphs is ~ 30%. It thus 
appears that the four dSphs as a group tend to overproduce the 
number of galaxies with /* distributions that significantly de- 
viate from the mean, suggesting that the observed cumulative 
metallicity functions are not identical. This can be interpreted 
as an indication of clustered star formation and that a signif- 
icant fraction of the observed stars (i.e., at [Fe/H] < —2.5) 
in these galaxies were once formed in a cluster environment 
similar to that of the present-day Galactic Disk. It is, how- 
ever, important to identify star clusters in multi-dimensional 



C-space, as with the cluster in Sextans, in order to firmly es- 
tablish the existence of a clustered star formation signature. 

We note that below [Fe/H] ~ —3 or so, the observations 
tend to show a smaller scatter than what is expected from the 
models. Although the observed sample of galaxies is admit- 
tedly very small, this could suggest that the formation process 
of star clusters in the lowest metallicity regime was different 
from that in the Milky Way today, given that most present-day 
clusters are formed in the designated mass range. A smaller 
dispersion between the cumulative functions is achieved, e.g. 
for a larger fraction of smaller clusters. However, due to the 
small number of stars and relatively poorly determined metal- 
licities, Starkenburg et al. (2010) noted that the exact shape 
of the observed /* should be considered with caution, in par- 
ticular below [Fe/H] = —3. Deeper observations and larger 
samples of galaxies are certainly required in order to success- 
fully probe the cluster formation process in this metallicity 
regime. 

5.1.2. Star formation in ultra-faint dwarf galaxies 

We also ran a series of simulations of ultra-faint dwarf galax- 
ies with the same setup as for the dSph-sized galaxies. We 
assume a fiducial luminosity of L ga \ = 1 x 10 4 Cq for 
all model galaxies. This translates to an initial mass of 
-Mgai,init ~ 4 x 10 4 Mq, given t hat t he galaxies only con- 
sist of old stars as in Sextans (cf. 33.21 . In general, a bigger 
scatter between the /* is therefore to be expected as a result 
of the much lower galaxy-to-cluster mass ratio relative to the 
dSphs. To allow a more direct comparison with the models 
of the dSphs, we assume that each ultra-faint MDF peaks at 
[Fe/H] ~ —2 (cf. Fig. |3}. This is probably at the high end 
for most ultra-faint dwarf galaxies in the Milky Way system, 
which appear to peak in the range —2.5 < [Fe/H] < —2 
(see Kirby et al. 2008). However, this should not affect the 
conclusions presented here. 

Several of the trends for the ultra-faints are the same as 
in the corresponding dSph case. As shown in Fig. [7] the 
ICMF with 7 = 1 (bottom left panel) produces cumulative 
metallicity functions that are substantially different from each 
other. Since the model galaxies have a very low stellar mass 
(-Mgal.init ~ 4 x 10 4 Mq), they contain at most a few clusters 
below [Fe/H] = —2.5 (i.e., for 7 = 1). In such an extreme 
case, the assumptions made for generating the /* are on the 
brink of breaking down, wherefore the exact shape of the /* 
should be regarded with some caution. Having said that, the 
cumulative metallicity function of Fornax (Fig. |6]l appears 
qualitatively similar to many of the model /* in the bottom 
left panel of Fig. [7] with a lack of stars below [Fc/H] ~ —3. 
This could suggest that Fornax was, to a high degree, pre- 
enriched in metals before forming the first generation of low- 
mass stars, possibly by a supermassive cluster. However, fur- 
ther studies of its metal-poor tail in combination with accurate 
modeling of the chemical evolution is called for in order to de- 
termine whether pre-enrichment occurred to a larger extent in 
Fornax than in other dwarf galaxies, and, if so, what would be 
the reason for this pre-enrichment. 

As for the dSph case, the ICMF with 7 = 2.5 (Fig. |7] mid- 
dle left panel) show no obvious signs of clustering, despite 
the low mass of the simulated dwarf galaxies. For 7 = 2 (top 
left panel), there is a clear signature of clumping with a ma- 
jor fraction of the /* considerably deviating from the mean 
(cf. Fig. [7] middle right panel). The same trend is present 
in the dSph case, only weaker. In contrast to the dSph case, 
however, the ICMF inferred from the simulations by Wise et 
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Figure 7. Predicted cumulative metallicity functions, /„ (normalized to unity at [Fc/H] = —2.5), for typical ultra-faint dwarf-sized galaxies (cf. Fig. |6). Each 
galaxy is assumed to have a present-day luminosity of £ ga i ~ 10 4 C@, implying an initial stellar mass of Af ga i init ~4x 10 4 A4@. Here, each /„ is sampled 
with 12 stars, reflecting the average number of stars in the survey of ultra-faint dwarf galaxies by Kirby et al. (2008, 2010, 2011b). The sand-shaded regions, 
indicating the typical Poisson noise due to the finite sampling, are changed accordingly. In the bottom right panel, 8 ultra-faint dwarf galaxies observed by Kirby 
et al. (2008, 2010, 2011b) are displayed. These are: CVn I (down-faced triangles), CVn II (circles), UMa I (squares), UMa II (up-faced triangles), Leo IV 
(diamonds), ComB (left-faced triangles), Here (hexagons), and Leo T (pentagons). Otherwise, symbols are the same as in Fig. [6] 



al. (2012b) appears to generate an increased number of /* 
that deviate from what is expected from the finite sampling of 
12 stars (Fig. [7] top right panel). Interestingly, the observed 
/* (Fig- El bottom right panel) of 8 ultra-faint dwarf galaxies 
(Kirby et al. 2008, 2010, 2011b) show no detected disper- 
sion that can readily be attributed to the formation of clusters. 
Only three out of eight /„ are at some point found outside 
the shaded areas, again taking into account the difference in 
the sampling of the Statistically, this should occur 1.6 out 
of eight times, if the cumulative metallicity functions were 
drawn from a single parent distribution (i.e., no clustered star 
formation). Hence, it appears that clumping is underproduced 
in the ultra-faints, also as compared to the case in which the 
ICMF of Wise et al. (2012b) is assumed. This is at odds with 
what naively would be expected. In a set of simulations, Wise 
et al. (2012a) investigated the impact of momentum transfer 
from ionizing radiation. The ICMF in these simulations was 



found to have a slope of 7 ~ 3.5 — 4 (J. Wise, priv. comm.). 
Such a steep slope would significantly reduce the number of 
more massive clusters, which are mainly responsible for the 
excess clumping in the /* (cf. Fig. [7] middle left panel), in 
accordance with observations. However, why the radiation 
pressure in that case would not inhibit the formation of mas- 
sive clusters in the dSphs remains to be understood. 

One explanation to the apparent overproduction of irregu- 
larities in the top panels of Fig. [7] as compared to observa- 
tions is that we have assumed too low stellar masses in the 
modeling of these galaxies. The ultra-faints could have been 
stripped off most of their stellar mass due to tidal interactions 
with the Milky Way. If so, the number of star clusters that 
constituted the original galaxies would have been higher, and 
each cluster would therefore have been sampled with fewer 
stars, leading to smoother However, from N-body sim- 
ulations of the evolution of dSphs driven by galactic tides, 
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in Figs.|6]and|7] 



Penarrubia et al. (2008) conclude that ultra-faints are unlikely 
to be rest products of tidally stripped, more massive dwarf 
galaxies. Tidal stripping predicts to steep a change in the 
velocity dispersion as compared to what is observed. This 
finding is in line with the existence of metallicity gradients 
in several of the dSphs (de Boer et al. 2011), which indicate 
that stripping due to tidal interactions at later stages could not 
have been overly severe. Furthermore, the observed Z — L re- 
lation for dwarf galaxies (e.g. Kirby et al. 2008), which states 
that the mean metallicity of the stellar population decreases 
with decreasing luminosity (i.e., stellar mass), argues against 
such a scenario. Another alternative is that our assumption 
on the chemical homogeneity of star clusters may be incor- 
rect. However, we have no a priori reason to believe that the 
physical processes responsible for the homogenization were 
significantly different in the early universe. Besides, the iden- 
tification of the three potential cluster stars in Sextans (Fig. [TJ 
argues for the formation of chemically homogeneous clusters, 
also at early times. 

We conclude that the lack of pronounced irregularities in 
the cumulative metallicity functions of the ultra-faint dwarf 
galaxies is a signature of reduced clustering. Evidently, the 
bulk of the metal-poor clusters that formed in ultra-faints had 
masses (i.e., the fraction of mass locked-up in low-mass stars) 
below the current detection limit. This seems to be in contrast 
to the dSphs, which, despite their significantly hig her (st ellar) 
masses, show some evidence of clustering (see £15.1. 11 1. We 
will further discuss the possible implications of this finding in 
the next section. 



5.2. Implications for galaxy formation 

Before we immerse into the intricacies of galaxy formation, 
let us have another look at the data of the dSphs and the ultra- 
faints. As an additional check of our findings in £15.11 we com- 
puted the mean /*,dSph and /^ultra-feint from the observa- 
tional data of the dSphs and ultra-faint galaxies, respectively 
(Fig. [8] lower panels). These averages have the advantage 
of being model independent. However, due to the relatively 
low sampling, they tend to follow the irregularities of the in- 
dividual galaxies. As a consequence, the chance is higher that 
a given irregularity of an /* is found inside the region de- 
fined by the scatter due to the finite sampling (shaded areas). 
The probability of detecting pronounced clumping is there- 
fore slightly decreased. Certainly, larger samples of galaxies 
would be desirable, especially a larger sample of dSphs. With 
this in mind, we note that two out of four dSphs are found out- 
side the shaded area in Fig. [8](l ower le ft panel). Although this 
is fewer than what is found in £15.1.11 it is more than what is 
expected from the finite sampling alone and still in agreement 
with the excess scatter predicted for the fiducial ICMF (upper 
left panel). To summarize, the fraction of dSphs with detected 
excess scatter is 50 — 100%, depending on how the boundaries 
are defined. Recall that the corres ponding prediction for the 
fiducial ICMF is 75% (see £15.1. 11 1. In contrast, only two out 
of eight ultra-faint dwarf galaxies (and one border-line case) 
are found outside the shaded area i n the lower right panel of 
Fig. [8] in line with the result in £15.1.21 Hence, given that 
cluster masses were distributed according to the present-day 
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ICMF also in the ultra-faints (Fig. [8] upper panels), the ex- 
pected scatter is noticeably larger than what is observed. Evi- 
dently, the fiducial ICMF overproduces the amount of scatter 
and does not seem to adequately describe the birth distribu- 
tion of star clusters in the ultra-faint dwarf galaxies, unlike in 
the dSphs. 

Figure [8] also reveals another distinction between the two 
dwarf galaxy populations. The average cumulative metallic- 
ity function of the ultra-faints (connected, green diamonds in 
the lower right panel) rises jibove the corresponding cumu- 
lative metallicity function, /*,Hak» observed in the Galactic 
Halo (Schorck et al. 2009; connected, black diamonds). In 
contrast, /*,dSph (connected, red diamonds in the lower left 
panel) closely follows /*,Halo- The difference is quite notice- 
able. The horizontal error bars on the dwarf galaxy data in- 
dicate the la uncertainty. At face value, it appears that the 
fraction of extremely metal-poor stars in the classical dSphs, 
on average, is lower than in the ultra-faint galaxies, while 
similar to that of the Galactic Halo. Note that if, in fact, 
/*,dSph = /*,Haio, all dSphs in the current sample will at 
some point be fou nd outside the shaded area as defined by 
/*,Haio (see i j5.1.U . 

Is the difference between /*, d s P h and /*, u itra-feint real? 
We note that Fe abundances of the stars in the dSphs are es- 
timated from the strength of the CaT, while stellar Fe abun- 
dances of the ultra-faints are determined directly from a large 
number of metallic lines, including Fe. A suppression of 
the metal-poor tail of the MDF may arise, e.g. if the CaT 
method is poorly calibrated, a trend in the [Ca/Fe] ratio with 
[Fe/H] exists that is not properly accounted for, or if there ex- 
ists a pristine, extremely metal-poor population of stars in the 
dSphs that is less centrally concentrated and which therefore 
is not fully sampled. The CaT method have recently been 
re-calibrated to hold down to [Fe/H] ~ —4 (Starkenburg et 
al. 2010). Furthermore, no trend in [Ca/Fe] is observed be- 
low [Fe/H] = —2.5, neither in the Galactic halo (Cayrel et 
al. 2004), nor in the dSphs (e.g. Koch et al. 2008; Cohen 
& Huang 2010; Tafelmeyer et al. 2010; Kirby et al. 2011a), 
although the present data in this metallicity regime are admit- 
tedly very sparse for the dSphs. Kirby et al. (2011b) have 
also observed MDFs for a number of dSphs, including Sex- 
tans, Sculptor, and Fornax. Examining their data, we see 
no obvious deviation from the observations by Starkenburg 
et al. (2010). Metallicity and age gradients could be a po- 
tential problem, where the oldest and most metal-poor pop- 
ulations usually are found in the outskirts of dwarf galaxies. 
Battaglia et al. (2011) also detect a metallicity gradient in 
Sextans. However, judging from their Fig. 8, it is seems to 
flatten out towards R = 2°. Besides, a less centrally con- 
centrated, pristine stellar population in the dSphs must have 
a peak metallicity below [Fe/H] ~ — 3 (and be accompanied 
by the absence of such a bias in the ultra-faints) in order to 
compensate for the lack of extremely metal-poor stars. This 
is not clearly the case, at least not for Sextans. Furthermore, 
as regards /* i Haio, Schorck et al. (2009) note that the rejection 
of stars with strong G-bands in the Hamburg/ESO data might 
lead to a bias, in particular against stars with [Fe/H] < —3.5. 
They conclude, however, that a possible effect on the inferred 
MDF would only be minor. At this point, we will take the cu- 
mulative metallicity functions of the dSphs, ultra-faints, and 
of the Galactic Halo as given. 

5.2.1. Different origin of the dSphs and the ultra-faints? 



Based on the observed difference in the degree of clumping 
in [Fe/H]_between dSphs and ultra-faints and on the offset 
between /*, d s P h and j*, u itra-faint, we argue that the two pop- 
ulations of galaxies have formed in different environments, in 
which the conditions for star formation were not the same. In 
their on cold dark matter (CDM) cosmological hydrodynam- 
ics simulations, Bovill & Ricotti (201 la, b; see also Bovill & 
Ricotti 2009; Ricotti 2010) identified an early population of 
faint dwarf galaxies that were able to form stars only prior to 
reionization (except a fraction of galaxies that may have ex- 
perienced late-time gas infall and subsequent star formation). 
These "pre-reionization fossils" are defined as having a max- 
imum circular velocity remaining below v c = 20 km s _1 at 
all times and a luminosity < 10 5 Cq. The authors argued that 
at least a fraction of the ultra-faints are consistent with being 
"pre-reionization fossils". Salvadori & Ferrara (2009) found 
in their hierarchical merger tree simulations that the simulated 
galaxies with similar luminosities and mean metallicities as 
the ultra-faints were formed within H2 -cooling minihaloes at 
z > 8.5, i.e., before reionization. This result is in accord with 
that of Bovill & Ricotti (2011a, b). Star formation in these 
minihaloes is supposed to be critically affected by the inef- 
ficient cooling, including a higher cosmic microwave back- 
ground (CMB) temperature, and be susceptible to feedback 
effects like H2 -dissociating Lyman-Werner radiation. The 
characteristic cluster mass could therefore have been much 
lower than in the present-day Galaxy (cf. Stacy et al. 2010; 
Clark et al. 2011), and still lower than in the dSphs, if this 
population of galaxies predominantly were assembled after 
reionization. 

The reason for the offset between j*,ds P h and /* :U itra-faint 
is not clear. Seemingly, the dSphs exhibit an effect similar 
to the classical "G-dwarf problem" in the Solar neighborhood 
(van den Bergh 1962). In our present stochastic chemical evo- 
lution model (see also Karlsson 2006), we make use of the 
Z?trans criterion for low-mass stars formation (Frebel et al. 
2007). This means that low-mass stars can only be formed in 
carbon and/or oxygen enhanced gas for which the criterion: 
D trans = log 10 (10l c / H ] + 0.3 x 10[°/ H 1) > -3.5 is ful- 
filled. This effectively suppresses the number of stars below 
[Fe/H] ~ —3.5 in the model, mimicking an IMF with a higher 
characteristic mass. Interestingly, the prediction is in close 
agreement with the data of the ultra-faints (Fig. [8] lower right 
panel). The suppression is, however, not enough to explain the 
relative lack of extremely metal-poor stars in the dSphs, why 
some kind of additional suppression must have occurred. Fol- 
lowing the argument above, an offset between the cumulative 
metallicity functions could be induced if the enrichment of 
carbon, relative to Fe, was lower in the dSphs than in the ultra- 
faints. The fraction of carbon-enhanced metal-poor (CEMP) 
stars (i.e., those with no s-process enhancements) may give us 
a hint whether different levels of C enrichment were a reality 
at early stages in the two types of galaxies. It is worth noting 
that several CEMP stars have recently been discovered in the 
ultra-faints (Norris et al. 2010a, b; Lai et al. 2011). An alter- 
native explanation to the offset is that the dark-matter haloes 
in which the more luminous dSphs reside, were accreting gas 
from the intergalactic medium that was already slightly pre- 
enriched in metals. This scenario is in line with the results 
of Salvadori & Ferrara (2009). A continuous infall of suffi- 
ciently metal-poor gas would further enhance the effect (cf. 
Kirby et al. 2011b). 

From simulations of primordial stars (Bromm et al. 2009; 
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Karlsson et al. 2011 and references therein), it is widely ac- 
cepted that the special conditions in the Universe at the end 
of the cosmological "Dark Ages" led to a top-heavy primor- 
dial/early IMF. If these conditions prevailed in the ultra-faints 
more so than in the dSphs, we would expect this to be reflected 
in the stellar chemical abundance data. A close inspection of 
the [Mg/Fe] ratio in Fig. [2]may, indeed, hint towards such a 
difference. The mean [Mg/Fe] ratio below [Fc/H] = —2.5 
for the compiled dSph data (dark-grey symbols), excluding 
the data of Aoki et al. (2009), is ([Mg/Fe]) = 0.30 ± 0.14. 
The corresponding ratio including all the stars of the ultra- 
faints (light-grey symbols) is ([Mg/Fe]) = 0.56 ± 0.27. The 
la deviation denote the total deviation, including both ob- 
servational uncertainties and cosmic scatter. This is a fairly 
noticeable difference. At face value, a higher mean a/Fe ra- 
tio in the ultra-faints could be suggestive of a top-heavy IMF 
favoring the formation of more massive stars and would be 
in line with our conclusion above that ultra-faints and dSphs 
are distinct. However, the mean [Mg/Fe] ratio for the ultra- 
faints is based on 9 stars only. It is important for a comparison 
like this to be internally consistent to minimize systematic un- 
certainties, to avoid observational biases, and to have control 
over contamination effects such as SNe la enrichment, galac- 
tic winds, etc. Before more data are available, this particular 
finding should therefore be regarded merely as tentative. 

5.2.2. Remarks on the origin of the Galactic Halo 

Based upon their data, Aoki et al. (2009) conclude that dwarf 
galaxies like Sextans could not be a major contributor to the 
Galactic stellar Halo. In contrast to the bulk of metal-poor 
Galactic Halo stars (as well as to dSph stars, see Fig. [2j, the 
majority of stars in their very metal-poor sample is not en- 
hanced in a elements as compared to Fe. However, in the 
light of the results presented above, we would like to add a 
small caveat to this conclusion. If we accept the scenario that 
stars could form in clusters already at [Fc/H] ~ —3, and that 
the three "blue" stars once belonged to such a cluster, the ex- 
istence of subsolar [Mg/Fe] stars below [Fe/H] = —2.5 in 
Sextans could naturally be explained as a result of chemical 
inhomogeneities in the star-forming gas (see Fig. |5), with- 
out invoking changes in global parameters like the IMF. For 
that reason alone, dwarf galaxies like Sextans should not au- 
tomatically be rejected as "progenitors" of the Galactic Halo; 
they may just be statistical outliers. Having that said, we note 
that the low [Mg/Fe] rati o of the potential cluster still needs 
to be explained (see fj43|. More data are needed, particularly 
around —3 < [Fe/H] < —2 in order to determine whether a 
low [a/Fe] ratio is a global property of the galaxy already at 
these metallicities. 

Current data from the Sloan Digital Sky Survey (SDSS) 
suggest that the stellar Galactic Halo is composed of (at 
least) two components; an inner Halo with a peak metallic - 
ity at [Fe/H] ~ — 1.6 and a more metal-poor, net retrograde- 
rotating outer Halo with a peak metallicity at [Fe/H] ~ — 2.2 
(Carollo et al. 2007, 2010; see also, e.g. Schuster et al. 2012 
for evidence of a multi-component Halo at [Fe/H] > — 1). 
Carollo et al. (2007) speculated that the outer Halo was 
formed through dissipationless, chaotic merging of small sub- 
systems, such as faint dwarf galaxies. Recently, Schonrich et 
al. (201 1) claimed to have identified significant biases in the 
distance estimates of the stars in the Carollo et al. (2007, 
2010) sample, in particular of the lowest metallicity stars. 
They conclude that, as a result of these distance overesti- 
mates, the evidence for the existence of an outer and inner 



Halo is weakened (see, however, Beers et al. 2012). If ultra- 
faints were a major contributor, exclusively to an outer stel- 
lar Halo, we expect the cumulative metallicity function of the 
outer Halo to differ from that of the inner Halo. Moreover, if 
our findings hold, the outer Halo may also exhibit less clump- 
ing in C-space than the inner Halo (Karlsson et al. 201 1). An 
extensive, multi-elemental abundance survey of the Galactic 
Halo would be able to disentangle potential components and 
unveil the progenitor systems from which the Halo was once 
assembled. 

6. FINAL REMARKS 

In the published abundance data for six very metal-poor stars 
in the Sextans dSph by Aoki et al. (2009), we have identi- 
fied three stars around [Fe/H] = —2.7, that potentially were 
once members of a now disrupted star cluster. If so, this is 
currently the most metal-poor cluster known to have formed 
in the Universe. We estimate the initial mass of the cluster to 
M,,i„it = 1.9±J;9 (1.6±oJ) x 10 5 M e , assuming a Salpeter 
(Kroupa) IMF. Aoki et al. (2009) infer a very low [Na/Fe] ra- 
tio for one of the three potential cluster stars. This may be sug- 
gestive of an Na-O anticorrelation present in the data, imply- 
ing that the relic cluster was a GC. One out of three stars with 
low [Na/Fe] is typical for GCs (Carretta et al. 2009) and the 
existence of a relic GC at [Fe/H] = —2.7 in Sextans would 
also resonate with the presence of metal-poor GCs in other 
dSphs (Letarte et al. 2006). We argue that the metal-poor 
environment of dwarf galaxies like dSphs provides an ideal 
base for searching for ancient clusters. We note that "bumps" 
or irregularities in the metal-poor tails of the MDFs of dSphs 
could be a signature of clusters. In the data of Kirby et al. 
(201 la), several MDFs exhibit "bumpy" tails, such as Leo II, 
Draco and Ursa Minor. Interestingly, the Ursa Minor dSph is 
known to possess a second concentration of stars roughly 14' 
from the center. Kleyna et al. (2003) claimed that these stars 
belong to a cold substructure, consistent with the remnants of 
a disrupted star cluster. Whether the stars in the "bump" at 
[Fe/H] ~ —3 (Kirby et al. 2011b) have the same origin as 
the stars in the cold substructure remains to be seen. In their 
sample of ten Ursa Minor stars, Cohen & Huang (2010) de- 
tect two stars around [Fe/H] = —3.1. The abundance data is, 
however, inconclusive regarding their origin. 

Inspired by the possible detection of a metal-poor star clus- 
ter in Sextans, we present a new way of interpreting cumu- 
lative metallicity functions of dwarf galaxies to study the star 
formation process at high redshifts and in the most metal-poor 
regimes. A comparison between models and available metal- 
licity data for dSphs and ultra-faint dwarf galaxies suggests 
that the formation of clusters in the ultra-faints differed from 
that of the dSphs. Together with a detected offset between 
the mean cumulative metallicity functions of the ultra-faints 
and the dSphs, and a putative offset in the ([Mg/Fe]) ratio, 
we speculate that these two populations of galaxies probably 
were formed in different environments and would therefore 
be distinct in origin. A possible explanation to the apparent 
absence of clumping in the ultra-faints below [Fe/H] = — 2.5 
is that these galaxies were formed predominantly before the 
local Universe was reionized, while the dSphs formed mostly 
after reionization. However, whether there is a dichotomy, as 
suggested by the current observational data, or a continuous 
distribution of galaxies exhibiting a mix of properties should 
be investigated further. Recently, Lai et al. (2011) observed 
25 stars in the Bootes I ultra-faint dwarf galaxy, whereof 12 
stars were found below [Fe/H] = —2.5. The cumulative 
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metallicity function inferred from their data falls right on top 
of the mean (see Fig. [8] lower right panel) as defined by the 
data of Kirby et al. (2008, 2010, 201 lb), with no features that 
can be attributed to clustering. This is entirely in line with our 
conclusion that the environment in the ultra-faints seemed to 
have inhibited the formation of massive clusters, at least in the 
regime [Fe/H] < —2.5. A search for clumping in the ultra- 
faints above [Fe/H] = —2.5 should be performed in order to 
determine the impact of the metallicity on the ICMF. Finally, 
we note that if more massive clusters were extremely rare in 
the ultra-faints, this could potentially have a significant im- 
pact on the escape fraction of ionizing radiation from these 
galaxies, thereby influencing on their capacity as sources of 
reionization. 

Our new work has important implications for near-field cos- 
mology, in particular, in identifying unique chemical signa- 
tures in the local universe that we can attribute to the first 
stellar generations. In principle, it is possible to probe and to 
reconstruct ancient star clusters that have long since dissolved 
into the dwarf population. The reconstructed systems can tell 
us about the onset of stellar clustering and the build-up of the 
chemical elements with cosmic time. 

We believe that reliable detections of stellar clustering 
within old stellar populations will be commonplace in an era 
of extremely large telescopes. At the present time, while sev- 
eral multi-object facilities are under discussion, the only con- 
crete proposal is the MANIFEST project that will provide 
high resolution spectroscopy over the 20' field of the Giant 
Magellan Telescope (Saunders et al. 2010). This field of view 
is well matched to the size of largest dwarf galaxies in the 
Galactic halo. Bland-Hawthorn et al. (2010a, their Figs. 7 
and 8) present simulations of C-space clustering in low mass 
dwarfs and show that even a modest multi-object spectrograph 
will provide a tremendous advance on the capabilities we have 
today. It will be possible to obtain high quality abundances for 
stars down to r = 22 mag. This will allow us to extend our 
reach and observe dwarf giants to beyond 100 kpc in Galac- 
tocentric radius. 

Many more dwarf spheroidals and ultra-faint dwarf galaxies 
are likely to be identified by future all sky survey telescopes, 
including the Large Survey Synoptic Telescope. With enough 
resolved clusters in individual dwarfs, or treated collectively, 
we can learn about the star formation process before, during 
and after the epoch of reionization. These new data will re- 
quire a more rigorous theoretical framework than exists at the 
present time. Towards this end, it will be necessary to resolve 
star clusters in cosmological simulations from the onset of the 
first stars through reionization to the present epoch. 
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